Disorder effect in low dimensional superconductors 
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The quasiparticle density of states (DOS), the energy gap, the superfluid density p s , and the localization effect in 
the s- and d-wave superconductors with non-magnetic impurity in two dimensions (2D) are studied numerically. 
For strong (unitary) scatters, we find that it is the range of the scattering potential rather than the symmetry 
of the superconducting pairing which is more important in explaining the impurity dependences of the specific 
heat and the superconducting transition temperature in Zn doped YBCO. The localization length is longer in the 
d-wave superconducting state than in the normal state, even in the vicinity of the Fermi energy. 
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In this paper we present numerical results for 
the density of states and other physical quanti- 
ties in superconducting states with strong poten- 
tial scattering from non-magnetic impurities in 
Zn doped YBCO. A comparison with experiments 
as well as other theoretical results is also given. 

The model we study is the BCS mean-field 
Hamiltonian with disorder which is defined by 



H =-t clc lla + J2( A rrC^cl +Ti + h.C.) 
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where < rr' > denotes nearest neighbors, /j is 
the chemical potential, and r imp are the positions 
of random impurities. We assume the random 
potential V(r - r imp ) = V S TtTimp + Vi8 <rYimp> . 
When Vi = 0, it is a (5-function potential, oth- 
erwise it is finite ranged [1]. A rT is the super- 
conducting order parameter. It should be deter- 
mined self-consistently from the relation A rr = 
J < c r -fC I+T i >, where J (assumed disorder inde- 
pendent) is the coupling constant. In the absence 
of disorder, A rr is translation invariant with a 
symmetry compatible with the point group of the 
lattice. In this paper only the on-site s-wave pair- 
ing symmetry, A rr = A6 Tt o, and the d-wave pair- 
ing symmetry, A rr = A(<$ Ti ±£ — 6 T ,±y), will be 
considered. 

In the presence of disorder, A rT is no longer 
translationally invariant. However, if the fluctu- 
ation of A rr in space is not too strong, the sym- 
metry of the gap function is preserved on a suf- 
ficiently large scale. To simplify the calculation, 
we approximate A rT in (1) by their average in 
space, A r = J2 r A IT /N (N the lattice size), so 
that only A T is determined self-consistently. The 



results obtained from this approximation can be 
taken as a qualitative guide to the properties of 
superconductors in the limit of very strong disor- 
der, but are accurate for dilute impurities. 

Fig. 1 shows the DOS [2,3] for the s- and d- 
wave superconductors with strong impurity po- 
tentials for various impurity concentrations x ob- 
tained using a recursion method. For both sym- 
metries the DOS near zero energy increases with 
x and becomes comparable with the average DOS 
at x~0.07 in the finite range model. This behav- 
ior agrees qualitatively with the variation of 7(0) 
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FIG. f . DOS as a function of energy for the s- and 
d-wave superconductors at half-filling. 



(linear specific heat coefficients) with Zn concen- 
trations in the Zn doped YBCO [4]. For the d- 
wave superconductor with a strong (5-function po- 
tential, a small DOS peak is observed near the 
zero energy in consistent with the resonant scat- 
tering theory of the d-wave superconductor in the 
unitary scattering limit [5]. As the impurity po- 
tential becomes small, this peak disappears, but 
the DOS at the zero energy is always found to be 
finite. We have not found evidence for the sin- 
gularity at the zero energy which was predicted 
recently by Nersesyan et al [6]. 

In Figure 2 we compare the energy gap and 
superfluid density p s as functions of x for the 
s- and d-wave superconductors at zero tempera- 
ture [3] . p s is obtained by directly calculating the 
current-current correlator in the ground state. A 
finite range impurity potential suppresses both A 
and p s more strongly than a (5-function potential. 
Comparing the result with the drop of the super- 
conducting transition temperature with x in the 
Zn doped YBCO [7], we find that our result for 
both the s- and d-wave superconductors with a 
finite range impurity potential can give a qual- 
itative fit with the experiment, while the result 
for the s- and d-wave superconductors with the 
(5-function potential decreases much too slowly to 
fit the experiment. This is in agreement with our 
estimate for the scattering phase shift from the 
resistivity measurements [1]. 
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FIG. 2. Normalized energy gap and superfluid den- 
sity vs x for the s- and d-wave superconductors at half 
filling on a 10x10 lattice at zero temperature. 



We have also examined the localization ef- 
fect in the d-wave superconductor by calculat- 
ing the localization length of the quasiparticle 
using the scaling method. We find that the lo- 
calization length in a superconducting state is 
far longer than that in a corresponding non- 
superconducting state for all energies, hence the 
localization effect is much weaker in a supercon- 
ducting state than that in a normal state, and 
is unlikely to influence the behavior of p s for ex- 
ample. The scaling function for the d-wave su- 
perconducting state is the same as for the cor- 
responding non-superconducting state consistent 
with a single parameter scaling theory. 
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